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Abstract 

The effects of interferon-beta (IFN-P), one of the key immunotherapies used in multiple sclerosis (MS), on peripheral blood 
leukocytes and T cells have been extensively studied. B cells are a less abundant leukocyte type, and accordingly less is 
known about the B cell-specific response to IFN-p. To identify gene expression changes and pathways induced by IFN-p in B 
cells, we studied the In vitro response of human Epstein Barr-transformed B cells (lymphoblast cell lines-LCLs), and validated 
our results in primary B cells. LCLs were derived from an IVIS patient repository. Whole genome expression analysis identified 
115 genes that were more than two-fold differentially up-regulated following IFN-p exposure, with over 50 previously 
unrecognized as IFN-p response genes. Pathways analysis demonstrated that IFN-p affected LCLs in a similar manner to 
other cell types by activating known IFN-p canonical pathways. Additionally, IFN-p increased the expression of innate 
immune response genes, while down-regulating many B cell receptor pathway genes and genes involved in adaptive 
immune responses. Novel response genes identified herein, NEXN, DDX60L, IGFBP4, and HAPLN3, B cell receptor pathway 
genes, CD79B and SYK, and lymphocyte activation genes, LAGS and IL27RA, were validated as IFN-p response genes in 
primary B cells. In this study new IFN-p response genes were identified in B cells, with possible implications to B cell-specific 
functions. The study's results emphasize the applicability of LCLs for studies of human B cell drug response. The usage of 
LCLs from patient-based repositories may facilitate future studies of drug response in IVIS and other immune-mediated 
disorders with a B cell component. 
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Introduction 

Aberrant B cell activity is an underlying pathological mecha- 
nism in a range of immune-related disorders, from diseases that 
are primarily antibody mediated, such as systemic lupus 
erythematosus, to diseases considered to be primarily T cell 
driven, such as multiple sclerosis (MS) [1,2]. Evidence for the role 
of B cells in MS includes the oligoclonal IgG bands present in the 
cerebrospinal fluid of most MS patients; the ectopic B cell follicles 
within germinal centers in meninges of MS patients; and the 
aberrant B cell cytokine response observed in MS patients [3] . The 
Epstein-Barr virus (EBV) that targets and infects B cells has been 
suggested as an MS risk factor, with recent supportive evidence 
including the increased disease severity in animals infected with 
the virus and the enrichment of EBV-infected B cells in MS lesions 
[4,5] . The pivotal role of B cell in MS was further demonstrated 
following reports of the therapeutic efficacy of the B cell-depleting 
anti CD20 drug in MS clinical trials [6,7]. Both antigen presenting 



functions of the B cell and antibody secretion contribute to disease 
pathogenesis [8,9]. 

Whereas much attention has been allocated in recent years to 
investigate B cell functions in MS pathogenesis, specific informa- 
tion regarding the effects of the MS immunomodulatory drugs on 
B cell activity is still lacking. This is also true with respect to 
interferon beta (IFN-P), which is among the key first-line therapies 
approved for relapsing MS [10]. Gene expression analyses and 
studies of cellular signaling pathways activated in response to IFN- 
P mosdy described the response in T cells, or peripheral blood 
leukocytes [11-13]. The IFN-P responses of small, yet important, 
cell populations such as B cells or monocytes have not received 
adequate attention. Indeed, accumulating evidence indicates 
immune cell-specific effects in gene expression are rife, including 
specifically for IFN-P, as shown by us and others for distinct genes 
and pathways in monocytes, in comparison to T cells [14-16]. 

Since B cells participate in both innate and adaptive immune 
functions, functions which are modulated by IFN-P [17] and 
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contxibute to MS pathogenesis, it is of interest to delineate the 
cellular pathways triggered by IFN-P in B cells of patients with 
MS. Previous studies described a role for IFN-P in B cell 
maturation [18], and in modulation of the TLR7 signaling 
pathway, a pathway involved in viral response and autoimmunity 
[19-21]. IFN-P was reported to increase transcription of the pro- 
survival genes encoding the B cell activating factor {TNFSF13B/ 
BAFF), STAT3, and STAT5 [15,22], but also to suppress release 
of stimulator)' cytokines such as IL-ip and IL23 [23]. In all, these 
studies focused on a few specific and known pathways, which 
suggest a complex manner in which IFN-P affects B cells. The 
discovery of other pathways, possibly B cell-specific, that are 
affected by IFN-P can benefit from a hypothesis-free approach, 
such as genome-wide gene expression analysis. 

The study of primary B cell functions is hampered by their low 
percentage among peripheral blood leukocytes and their short 
lifespan in culture. However, human B cells can be transformed by 
the EBV to create lymphoblast cell lines (LCLs) that retain many B 
cell properties, including typical cell surface markers and 
immunoglobulin secretion [24], as well as phenotypic character- 
istics of the donor B cells [25]. The usage of LCLs carries 
important technical advantages since they are much more easily 
maintained compared to primary B cell cultures and the in vitro 
set up allows repeated and parallel experiments on the same 
samples. Moreover, recent studies have confirmed that LCLs 
preserve the inter-individual gene expression variability of primary 
B cells, including heritable patterns governed by genetic 
background [26,27]. Despite the inherent limitations of research 
using transformed cell lines (reviewed in Welsh et al., 2009), it is 
increasingly recognized that LCLs derived from patient popula- 
tions can be useful for studies of the effects of genetic variation on 
cellular function as related to disease and drug response [28-32]. 
Studies probing the regulatory- functions of MS associated genomic 
regions have used genetic data from a LCL to demonstrate an 
overlap with B cell-specific transcriptionally active genomic 
regions [31,32]. Specifically, LCLs may serve as a model to study 
genetic variance related to the drug effects with specific reference 
to B cell activity. 

In this study, we used LCLs from people with MS to 
characterize the gene expression patterns of the response to in 
vitro IFN-P exposure. Whole genome analysis highlighted genes 
previously unknown to be affected by IFN-P. The results from 
analysis of LCLs, identifying new IFN-P response genes, were 
replicated in primary B cells, corroborating the relevance of the 
LCL system as a model for studying drug responses in B cells. 

Materials and Methods 

Study participants 

This study was appr()\'ed by the Carmel Medical Center 
Helsinki Committee and the Israeli Ministry of Health National 
Helsinki Committee for Genetics Studies. All participants provid- 
ed a written informed consent. Participants (>18 years) were 
recruited at the MS center at Carmel Medical Center, Haifa. 
Participants included were patients with clinically definite or 
laboratory supported MS diagnosis according to Poser and 
MacDonald criteria, [33,34] and healthy individuals as controls. 
Exclusion criteria for controls were presence of MS in family 
members up to third degree, or presence of any autoimmune or 
chronic inflammatory condition. Detailed demographic and 
clinical data were obtained from all subjects. Ethnicity was 
determined by the participant's self-report. MS disease type 
(relapsing-remitting, secondary-progressive, relapsing-progressive 
or primary-progressive) was recorded at date of phlebotomy. 



Blood samples were obtained from all participants, and peripheral 
blood mononuclear cells (PBMC) were purified using FicoU 
gradient (NovaMed). 

Lymphoblastold cell lines generation and culture 
conditions 

Frozen PBMC samples from participants were EBV trans- 
formed at the European Collection of Cell Cultures, England 
(http://www.phe-culturecollections.org.uk/3083.aspx), and the 
National Laboratory for the Genetics of Israeli Populations at 
Tel Aviv University, Israel (http://nlgip.tau.ac.il). The LCLs were 
cultured in RPMI-1640 supplemented with 10% FBS, 2 mM L- 
glutamine, penicillin (100 U/ml), streptomycin (100 |J,g/ml), and 
nystatin (12.5 U/ml, Biological Industries). LCLs were assayed 
within one month of thawing for the different experimental assays 
to minimize the passage numbers and avoid immortalization and 
other changes associated with long term growth in culture [25,35]. 

Flow cytometry 

LCLs were stained by multicolor flow cytometry with FITC 
anti-human CD 19; APC-CY7 anti-human CD27; PerCP anti- 
human CD38; and APC anti-human CD24 (Biolegend). BD 
CompBeads (Becton and Dickinson) were used for compensation 
according to the manufacturer's instructions. Unstained cells were 
used to exclude background fluorescence and isotype controls to 
determine antibody specificity. LCL viabilit}' was monitored using 
7-amino-actinomycin D (cBiosciences, USA). Data reading and 
acquisition was performed using the Flow Cytometer Cyan ADP 
(Beckman Coulter). Flowjo Software (TreeStar) was used for data 
analysis. 

Primary B cell Isolation 

B cells were isolated from PBMCS using the negative selection 
Human B cell isolation kit (Miltenyi), or EasySep Human B cell 
enrichment kit (StemceU) according to the manufacturers' 
instructions. All B cells preparations had at least 80% purity 
levels by CD 19 flow cytometry. 

IFN-p incubation assay 

Cells were incubated with 100 units/ml recombinant human 
IFN-Pla (InterferonSource) for 4, 16, or 48 hours. Cytokine 
concentrations used were within the range of reported serum IFN- 
P levels following drug injection in MS patients, and for the 
extended incubation times IFN-P was added to medium after 
24 hours to account for cytokine half-life [36,37]. Untreated cells 
were cultured in parallel and under identical culture conditions for 
comparison. LCLs and primary leukocyte IFN-P assays were 
performed at a density of 1.4 million cells/ml. 

RNA extraction 

RNA from LCL cultures was extracted using the Roche RNA 
isolation kit. RNA from primary B cells was extracted using the 
Ambion RNAqueous® Micro Kit. The concentration of RNA was 
measured by NanoDrop ND-100 Spectrophotometer (NanoDrop 
Technologies). RNA samples processed for the microarray chip 
assays and all RNA samples derived from primary B cells were 
checked for quahty using Experion (Bio-Rad). 

Microarray chip gene expression analysis 

A total of 16 pairs of RNA samples from LCLs (IFN-P treated 
for 4 hours or untreated) were hybridized to lUumina BeadChip 
HumanHT-12 v4, according to the manufacturer's protocols. Raw 
data was processed by GenomeStudio software (Illumina) for 
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quality control. Data were log2 transformed, and batch normal- 
ization was applied. Data were filtered for background noise and 
variance (variance S2° '). 

The gene expression dataset has been deposited in GEO with 
access no. GSE58240. 

Real Time RTPCR assays 

cDNA synthesis was performed using M-MLV reverse tran- 
scriptase (Promega) with random hexamer primers (Biological 
Industries). Specific PGR primers for selected genes were designed 
by us or obtained from published primer databases (PrimerBank, 
http://pga.mgh.harvard.edu/primerbank/ or qPrimerDepot, 
http://primerdepot.nci.mh.gov), and were obtained from Sigma 
or Integrated DNA Technology, IDT (Table SI). RT-PCR 
analyses were performed using the FastStart Universal SYBR 
Green Master-Rox (Roche) on the AB 7300 sequence detection 
system in duplicates. The ubiquitin-conjugating enzyme E2 D2 
{UBE2D2) gene was used as a reference gene as in previous and 
the present studies its levels were unaffected by IFN-P [,38]. 
Relative quantification of mRNA expression was calculated using 
the comparative GT method and is shown as fold change of 
expression (2"^^'^t) [39]. 

Statistical Analysis 

Statistical analysis of BeadChip gene expression array data was 
performed using JMP® Genomics version 6.0 software (SAS 
Institute Inc.). A total of 3035 probes out of the 47323 probes 
included in the BeadChip passed filtrations and were included in 
the statistical analysis. Principal component analysis and variance 
component analysis were performed as implemented in JMP 
Genomics. One-way ANOVA was performed to identify differ- 
entially expressed genes between IFN-|3 treated or untreated LCL 
samples, with correction for multiple testing by the False Discovery 
Rate (FDR) method [40]. Differentially expressed genes were 
defined as transcripts that had a log2(fold change of expression) > 1 
and an FDR corrected P-valuc<0.(),5. To select for novel genes 
that had not been previously recognized as IFN-P response genes, 
the differentially expressed gene hst was tested against the 
Interferome v2.01 database, using the parameters of log2(fold 
change)>l, and for selected down regulated genes also log2(fold 
change)< — 0.4 [41]. The IFN-P differentially expressed gene list 
was compared with the list of differentially expressed genes 
between LCLs and their counterpart B cells based on the data of 
Caliskan et al., with log2(fold change)>l and FDR <().()5 [26]. 

Data from RT-PCR experiments and flow cytometry analyses 
of cell subset distributions were analyzed using SPSS version 18.0 
(IBM), P-values £0.05 were considered significant, and Bonferroni 
correction was apphed as relevant. To compare between data from 
LCLs groups or with primary B cells and PBMCs, Mann Whitney 
U test was used. For studying the IFN-P effect by comparing 
untreated samples to treated samples, the Wilcoxon signed rank 
test was applied. 

Pathway Analysis 

Functional networks and canonical pathways affected by IFN-P 
treatment were explored with Ingenuity Pathways Analysis 
software (Ingenuity Systems). The data set containing the 
differentially expressed genes, their corresponding expression 
values, and P-values, was uploaded into Ingenuity software, and 
analysis was performed with a cutoff of P-value = 0.05 and 
log2(fold change) = 0.4 to obtain the top canonical pathways and 
functions. Top functional networks significantly activated (Z 
score>2) were further explored, as were specific upstream and 



downstream connections for genes of interest for their B cell- 
related function. 

Results 

B cell-subsets composition across donor groups is similar 

Both naive B cells and memory B cells express the EBV receptor 
CD21 and hence are the major targets of the viral transformation 
[42]. Thus, overall variance in gene expression across LCL 
samples may be influenced by differences in B cell-subset 
composition between samples. To address this issue B cell surface 
markers CD 19 and GD27 were analyzed in LCLs from healthy 
controls and MS donors, either treated with IFN-P or not at the 
time of phlebotomy, to evaluate B cell subset phenotype 
distributions [43] (Fig. 1). No differences were observed among 
the LCL groups in the major cell subset distributions representing 
naive B cell-like and memor)' B cell-like LCLs. Analysis of variance 
for B cell subset distributions revealed low variance within 
independent tests of a specific sample, but a significant variance 
between LCLs from different donors (ANOVA, P<0.001). These 
results indicate that the relative proportions of naive B cell-like and 
memory B cell-like cell subsets remain stable in culture for each 
sample, at least within the short-term culture period used in the 
current study. Moreover, we did not detect any segregation in cell- 
surface marker distribution that was related to the donor's disease 
status or therapy status at the time of phlebotomy. 

Gene expression analysis of IFN-pS response 

Whole genome gene expression analysis was performed for the 
in vitro response of 16 LCL samples to four hours IFN-P exposure. 
LCLs were obtained from donors affected with MS [63% females. 



100% n 




RRMS RRMS SPMS HC 
(untreated) (treated) 



Figure 1. Flow cytometry indicates similarity in proportions of 
cells with a na'i've B cell or memory B cell marker phenotype 
among LCLs from healthy donors and patients with MS. LCLs 
were typed for memory and naive B cell cell-surface phenotype by flow 
cytometry, by CD19 and CD27 immunoreactivlty. The percentage of 
cells for each subset Is shown for LCLs from donors with relapsing 
remitting MS that did not receive any therapy prior to the time of 
phlebotomy (RRMS: untreated, n = 12), donors with relapsing remitting 
MS under IFN-(3 therapy (RRMS: treated, n = 12), donors with secondary 
progressive MS (SPMS, n = 9), and donors that were healthy controls 
(HC, n = 5). Black bars: CD19+CD27- naive B cells; white bars: CD19+ 
CD27+ memory B cells. No differences were observed between all 
comparison groups (Mann Whitney test). 
dol:l 0.1 371/journal.pone.01 02331 .gOOl 
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Age median 38.5 yrs (range 20-59 yrs), of Jewish (50%) and Arab 
descent (50%), 81% with relapsing remitting MS, and 75% treated 
with IFN-P, at the time of blood collection]. Variance component 
analysis detected variance among donors and IFN-P response as 
major contributors to the overall variance in gene expression 
(Fig. 2). In addition, variance related to the donor state at time of 
phlebotomy, such as age, treatment status (whether donor was 
receiving IFN-P therapy or not), or disease subtype (relapsing 
remitting MS, secondary progressive MS or progressive relapsing 
MS) contributed each more than 5% to the overall variance. 

At our stringent cutoff, twofold change and above, at an 
adjusted P-value<0.05, 115 up-regulated genes were found 
following IFN-P treatment, but no down-regulated genes 
(Fig. 3A, Table S2 and Table S3). Cross-referencing the gene 
list with the Interferome V2.01 database [41] highhghted 55 genes 
that had not yet been reported as IFN-P response genes (Fig. 3B). 
Of these novel IFN-P response genes, 13 had been reported by 
Caliskan and colleagues as genes that are up-regulated in LCLs 
versus primary B cells and thus may be related to the viral 
transformation [26]. The remaining 42 genes may be considered 
as novel IFN-P response genes that are likely to be expressed in 
primary B cells as well, and from which four were selected for 
RTPCR validation as described next (3.3). 

Validation of novel IFN-|3 response genes in LCL and in 
primary B cells 

Four novel IFN-P response genes, NEXN, HAPLN3, DDX60L, 
and IGFBP4, were selected for additional validation based on 
their robust fold change values on one hand, but lack of functional 




0 0.05 0.1 0.15 02 025 0.3 0.35 



Weighted Average Proportion 

Figure 2. Proportion of variance in gene expression is 
explained by a major IFN-p response component and donor- 
specific components. Variance component analysis was modeled 
with the following components: Donor - representing variance 
contributed by unspecified differences between donors; IFN-[S exposure 
{in vitro IFN-p treated LCL sample versus the untreated sample); Age 
(stratified to under 40 years old or above); Gender; MS subtype 
(relapsing remitting, relapsing progressive, secondary progressive). 
Treatment status (IFN-p treatment naive donors, or donors treated for 
at least 1 year at the time of sample collection); and Ethnicity (Jewish/ 
Arab). The Residual component models all the variability that cannot be 
attributed to any of the explicit variance components. 
doi:1 0.1 371/journal.pone.0102331 .g002 



information with respect to B cell specific functions or connection 
to IFN-P related pathways. NEXN encodes nexilin, an F-actin 
binding protein proposed to function as a linker protein for the 
cytoskeleton and in focal adhesion junctions [44]. HAPLN3 
belongs to the hyaluronan and proteoglycan link protein family, 
and encodes a ubiquitous extracellular matrix protein [45]. The 
function of DDX60L can only be postulated from its homology to 
DDX60, an RNA helicase involved in the IFN antiviral response 
[46,47]. The insulin-like growth factor binding protein 4, IGFBP4, 
is a modulator of the activity of the insulin growth factors [48] . 

Using real time RTPCR these four genes were validated as IFN- 
P response genes, first in LCLs from MS patients and healthy 
donors, followed by validation in primary B cells and PBMCs 
(Fig. 4). We did not detect differences in the response to IFN-P 
between LCLs from MS and healthy controls for NEXN, 
DDX60L and IGFBP4, but for HAPLN3 a shght increase in 
expression was found in control samples (P = 0.048) (Fig. 4A). The 
expression of all four genes was significandy up-regulated (P^O.Ol) 
following IFN-P exposure in B cells, although to a different extent 
than in LCLs (P<0.03, Fig. 4B). An IFN-P transcription response 
for the four genes was also observed in PBMCs; however, for 
DDX60L and HAPLN3 the response was of a significandy lower 
magnitude than in B cells (P = 0.001), emphasizing the cell specific 
nature of this cytokine's response. 

To study the temporal dynamics of IFN-P, we assessed its effect 
on transcription of the new IFN-P response genes at longer 
exposure times up to 48 hours. All four genes selected for 
validation appeared to have a transient response to IFN-P, 
suggestive of early response genes, as at 16 hours of exposure their 
expression returned to baseline levels (Fig. 5). 

IFN-p induces activation of canonical IFN pathways and 
affects other immune pathways in LCLs 

Pathways analysis is especially useful for detection of coordi- 
nated expression changes of functionally related genes that 
individually do not pass the significance threshold. Ingenuity 
pathways analysis highlighted the known IFN response pathways, 
including the 'IFN signaling' pathway (P = 2.6E-13) and 'activa- 
tion of IRF by cytosolic pattern recognition receptors' (P = 1 .3E-6). 
In addition, the 'antimicrobial response' function is activated 
following IFN-P exposure (Z score = 2.8), shown in Figure 6, with 
the addition of nodes known to participate in the B cell innate 
response. The overall increased expression observed in this 
pathway is consistent with an IFN-P-induced activation of innate 
pathways in the LCLs. Thus, with respect to IFN-P canonical 
pathways, LCLs respond to IFN-P in a manner that is overall 
similar to other cell types. 

The top physiological systems affected were 'Hematological 
System Development and Function', and 'humoral immune 
response', both are systems that include functions of B cells (for 
both systems P<5.5E-03). Within the 'Hematological System 
Development and Function' category, the 'activation of lympho- 
cytes' network was increased (Z scores >2), including classic B cell 
genes like CD79B and TNFSF13B/BAEF, and genes tiiat are 
mainly known as T cell response regulators, like LAG3 and 
IE27RA [49-52] (Fig. 7A). The prototype B cell molecule CD79B 
that is part of the B cell receptor (BCR) complex and is expressed 
from early B cell developmental stages [53] was significandy down- 
regulated in LCLs (adjusted P = 0.01), albeit at a small fold change 
(0.7). Interestingly, many of the direct interacting partners of 
CD79B, including upstream transcription factors were down- 
regulated as well (Fig. 7B). Moreover, other components of the 
BCR pathway also decreased in expression following IFN-P 
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Figure 3. The LCL gene expression response to IFN-p. A: Volcano plot for the IFN-|3 gene expression response in LCLs. Log2(fold change) of 
expression levels following IFN-p exposure and P-values (one-way ANOVA) are depicted for each transcript probe. The horizontal dashed line 
represents the threshold for the adjusted P-value of 0.05. Vertical dashed line indicates threshold of log2(fold change)>1 and <-1. Differentially 
expressed genes (in black) are genes with Iog2(fold change) Sil and adjusted P-value <0.05, novel IFN-p response genes Included in RTPCR 
validation analyses are marked. For pathways analyses a lower threshold was employed (dotted vertical line) with log2(fold change)2:0.4 (dark grey 
dots). B: Venn diagram describing the proportion of differentially expressed genes in response to IFN-p [log2(fold change)2:1, adjusted P-value<0.05, 
bold outlined ellipse] that are novel response genes (grey shade), and a comparison to the published differentially expressed genes up-regulated in 
LCLs in comparison to primary B cells, based on the data of Caliskan and colleagues [log2(fold change)S:l, adjusted P-value<0.05][26]. Definition of 
novel IFN-p response genes is based on the Interferome V2.01 database (accessed 19.8.13) search for the 115 IFN-p response differentially expressed 
genes. Numbers indicate the number of genes within each subgroup. 
doi:10.1371/journal.pone.0102331.g003 



exposure, although the overall statistical significance did not reach 
the threshold (Fig. 7C). 

Several of the genes included in the lymphocyte related 
networks highlighted by Ingenuity pathways analysis are currently 
not recognized as IFN-P response genes by the Interferome 
database, including CD79B, SYK, IL27RA, LAG3, and KLF2. 
Quantitative RTPCR on RNA from primary B cells and PBMCs 
exposed to IFN-P confirmed these genes as IFN-P response genes 
(Fig. 7D). We note that for LAG3, up-regulation of expression in B 
cells was threefold higher in comparison to PBMCs, emphasizing 
the cell specific nature of the IFN-P effect. TNFSF13B/BAFF, 
which is not yet included in the Interferome database as an IFN-P 
response gene, was also up-regulated in B cells, in agreement with 
previous reports [15,22]. Thus, immune-related novel IFN-P 
response genes identified in LCLs were confirmed in primary B 
cells as IFN-P response genes. 

Discussion 

In this study we report that the cytokine and immunomodu- 
latory drug IFN-P activates prototype IFN canonical pathways in 
LCLs, including the antimicrobial innate immune response. Our 
novel findings comprise the decreased expression of several BCR 
pathway genes, including the immunoglobulin P component of the 
BCR CD79B together with its upstream regulators and down- 
stream interacting proteins following exposure to IFN-P. Thus, the 



LCL response to a short IFN-P exposure overall is similar to the 
known lymphocyte response, but it also includes a response 
component of t)'pical B cell genes and pathways. The effects of 
IFN-P on the BCR genes suggest a reduction of the capacity of the 
B cell to respond to antigen mediated signals, funneling its activity 
towards immediate innate system type of response. These results 
provide support to the applicability of LCLs as a model system to 
study the cellular effects of drug response, specifically with respect 
to B cell related functions. 

We identified in this study over 50 genes that were not 
recognized previously as IFN-P response genes using LCLs. The 
majority of these novel response genes have similar expression 
levels in LCLs and primary B cells [26] , and are therefore likely to 
be IFN-P response genes in B cells as well, as we were able to 
demonstrate for the genes NEXN, HAPLN3, DDX60L, and 
IGFBP4. Two of these genes, HAPLN3 a.nd DDX60L, displayed 
a large difference in the drug response between B cells and 
PBMCs, which may explain why they had not been previously 
identified as IFN-P response genes in studies using T cells or 
PBMCs. All four genes appear to be early response genes for IFN- 
P, and their transcription is under tight regulation, as indicated by 
the return to baseline levels after 16 hours of cytokine exposure. 
Presently, there is lack of functional data (as also indicated by lack 
of Ingenuity pathway connections) that can explain the contribu- 
tion of the activity of these four genes to IFN-P-induced B cell 
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Figure 4. Validation of tlie gene expression response for four 
novel IFN-P response genes by real time RT-PCit. The change in 
expression levels (2-'^'^'>) following four hours of IFN-(3 exposure is 
shown for four genes that are novel IFN-p response genes. The 
horizontal bars represent the median values. A: Comparison of the IFN- 
P response in LCLs from healthy controls (HC, empty symbols, n = 7) and 
patients with multiple sclerosis (MS, full symbols, n = 20). Fold changes 
for IFN-p response were significant for both groups (P-value<0.007, 
Wilcoxon signed rank test). No difference was observed in the IFN-p 
response between controls and MS LCLs, except for HAPLN3, for which 
the difference was of a small magnitude (*P-value<0.05, Mann Whitney 
test). B: The IFN-p response in primary B cells (n = 9) and PBMC (n = 8) 
from healthy controls and MS patients. All genes displayed a significant 
response to IFN-p for both B cells and PBMC (P-value<0.001, Wilcoxon 
signed rank test). Significant differences between the response in B cells 



and PBMCs is indicated by** P-value = 0.0003 (Mann Whitney test). The 
response magnitude of B cells and LCLs differed for all genes, however, 
these differences were small (P-value<0.03, Mann Whitney test). 
doi:1 0.1 371/journal.pone.01 02331. g004 

activity. NEXN and HAPLN3 may affect cell motility as 
cytoskeleton and extracellular matrix proteins, however, additional 
studies are required to assess their contribution to immune cell 
functional modulation by IFN-P. In the context of possible modes 
of action of IFN-P as an immunomodulatory therapy in MS, it is 
worth noting that IGFBPs have been suggested to Ik; involved in 
regulation of remyelination and increased levels of IGFBP4 in MS 
lesions were reported [54]. Elevated levels of IGFBP4 were 
reported to inhibit proliferation and induce apoptosis of mice 
lymphocytes in lymphoid tissues [48,55]. 

Some of the novel IFN-P response genes identified in this study, 
including KLF2, LAG3 and IL27RA, are known as immune 
function genes; however they had not been studied in the context 
of IFN activity and the effect on B cells. KLF2, the Kriippel-like 
factor 2, is a transcription factor involved in B cell homeostasis, 
regulating the development and balance between the marginal 
zone and follicular B cell pools [56]. KLF2 transcriptionally 
regulates the expression of several cell surface receptors that are 
involved in cell motility, however it appears to have a differential 
effect between B and T cells [56] . In B cells, KLF2 was shown to 
affect plasma cell homing to the bone marrow, mainly via 
transcription of integrin P7. Thus, down-regulation of KLF2 by 
IFN-P in B cells may decrease the migration capacity of antibody 
secreting B cells, affecting the adaptive response arm of the B cell 
function. LAGS is a CD4-like protein suggested to function as an 
inhibitory receptor to T cell activation and to regulate T cell 
homeostasis [50,57]. LAG3 is suggested to exert its inhibitory 
actions on T cells by competitive binding to MHC class II. The 
large increase in LAGS RNA levels mediated by IFN-P specifically 
in B cells, as shown herein, points to an inhibitory component on B 
cell activation that may be among the immunomodulatory effects 
of IFN-P. That the increase of LAGS levels by IFN-P may serve to 
attenuate cell activation in both B and T cells is in line with the 
modulatory role attributed to IFN-P in MS, the augmentation of 
autoimmunity by LAGS deficiency [58], and the increased 
activation of B cells occurring in MS [59]. IL27RA, a subunit of 
the receptor for IL27, is known to be required for maturation of B 
cells in the germinal centers through its signaling within the 
follicular T cells [51], it can induce IgG class switching in B cells 
[60], and is nxjuired for generation of high affmity anti-viral 
antibodies [52]. Thus, we suggest that down-regulation of IL27RA 
following IFN-P treatment may indicate the dampening of 
activation of B cell maturation, congTU(;nt with the effect of 
IFN-P on LAGS, decreased function of plasma cells due to KLF2 
down-regulation, and the overall decreased expression of the B cell 
receptor pathway genes. 

One of the important properties of the LCL as a research tool is 
that it can be used to study functional effects of genomic variation 
on expression, since LCLs recapitulate the donor variance in gene 
expression [26], although there is also contradictory evidence to 
this point [61,62]. This property has been recentiy exploited in 
association studies of gene expression and disease, and with drug 
response [6S-67]. In fact, Bolotin and colleagues demonstrate in 
their study the overall clustering of gene expression patterns of B 
cells and LCLs by donor [63]. Although our study was not 
designed to assess the effect of the donor's phenotype on the IFN-P 
response, our findings suggest that the gene expression variance 
among LCLs reflects, in addition to genetic variance, the donor's 
status at the time of phlebotomy, for characteristics such as age, 
disease state, and therapy status. The B cell-like subset composi- 
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Figure 5. NEXN, DDX60L, IGFBP4 and HAPLN3 respond transiently to IFN-p.The expression levels of NEXN, DDX60L, IGFBP4 and HAPLN3 in LCLs 
were analysed by real time RT-PCR analysis after 4 or 16 hours of exposure to 100 u/ml IFN-)3. Data is presented as a box plot of the 25"^ and 75"^ 
percentiles around the median values. 
doi:1 0.1 371 /journal.pone.01 02331 .g005 



tion of LCLs appears to be a sample-specific and reproducible 
characteristic that, similarly to gene expression patterns, may be 
related to the donor and his state at the time of phlebotomy. 
Studies employing animal models have been very instructive in 
deciphering immunological pathways, however, they have also 
revealed important differences in immune functions between mice 
and men, including B cell development and activity [68,69]. 
Moreover, IFN type 1 is species-specific [70,71], emphasizing the 



need to study the immunomodulatory properties of IFN- P activity 
in the context of human cells. Thus, analysis of drug response in 
experimental systems of animal models may fail to identify human- 
specific activity, and LCLs offer the advantage of an in vitro 
system that is based on human cells. Moreover, the study of 
complex drug interactions, including evaluation of drug mixtures, 
and comparative studies of drugs, can benefit from the in vitro 
experimental setting that allows controlled timing of drug 
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Figure 6. IFN-p activates the antimicrobial pathway and related innate immunity genes in LCLs. The pathway depicted includes the 27 
genes of the of the 'antimicrobial response' pathway that showed increased activation following IFN-p exposure of LCLs (activation score Z = 2.779) 
by Ingenuity analysis. In addition, genes outlined in blue were manually added to the pathway based on their known participation in the innate 
response in B cells or other leukocytes, and include TLR7 and TLR9 [20] [74], TRIM21 and TRIMS [75-77]. Only direct connections are featured, and 
subcellular location of gene products is as suggested by Ingenuity knowledge base. In green- down-regulated genes; in red- up-regulated genes, 
color intensity is proportional to the fold change in expression levels following IFN-p exposure. 
doi:1 0.1 371/journal.pone.01 02331. g006 
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Figure 7. Networks of B cell-related genes are affected by IFN-p in LCLs. Network of genes included in the 'Activation of Lymphocytes' 
function highlighted by Ingenuity analysis (Z = 2.33). Arrows points to IFN-p response genes that are not included in the Interferome database, and 
were included in a follow up validation analysis by quantitative RTPCR shown in D. In green- down-regulated genes; in red- up-regulated genes, color 
intensity is proportional to the fold change in expression levels following IFN-p exposure. A; An interaction network of CD79B, the immunoglobulin-p 
subunit of the B cell receptor signaling complex, with its regulators and interacting genes. The interacting nodes were limited to genes that were 
included in the filtered gene expression dataset, and were identified as directly connected to CD79B based on Ingenuity pathways analysis database. 
Notably, all the genes included in this network were found to be down-regulated by IFN-p (green) to some extent. Fold change and P-values are 
noted below each node. B: Decreased expression levels were observed for many of the genes included in the canonical B cell receptor signaling 
pathway in response to IFN-p. C: Quantitative reverse transcriptase RTPCR analysis of IFN-p response in primary B cells and PBMCs for genes marked 
in arrows in A and B. Genes were selected for validation because they were not identified by the Interferome database as IFN-p response genes, and 
are related to B cell function (CD79B, SYK, TNFSF13B/BAFF), or participate in the antimicrobial pathway depicted in Fig. 6 {IL27RA), or are related to 
activation of lymphocytes {LAG3, KLF2). * P-value = 0.007 (iVIann Whitney test). 
doi:1 0.1 371 /journal.pone.01 02331 .g007 
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application and dosage on one hand, while allowing consideration 
of die genetic and gene expression variance within the population, 
as expressed in the LCLs on the other hand. 

There are several known and inherent caveats in the in vitro 
system of the LCLs, which call for awareness and caution in data 
interpretation: First, the effect of the viral transformation on the 
gene expression patterns, and that of the viral programmed 
transcription on cytokine and drug induced pathways [26,63]. 
This is particularly pertinent when the study involves an antiviral 
cytokine or drug, like IFN-P. Secondly, the isolated and 
homogeneous cell culture m vitro lacks many of the functional 
interactions present in the in vivo multi-organ - multi-cellular 
environment, with regulatory feed-forward and feed-back effects 
emanating from the various cell types and their interactions. Thus, 
the LCL is an experimental system of lower complexity, with the 
advantages and disadvantages associated with a reductionist 
approach. Finally, time course of drug effects in the multicellular 
organism versus the cell culture is of an entirely different 
magnitude. The long term elfects of many drugs are cumulative 
and can differ significandy from the short term effects, as recentiy 
demonstrated for IFN-P in MS [72]. With reference to the present 
study, our study design led to identification of many genes which 
displayed transient changes in expression levels and returned to 
baseline levels at longer IFN-P exposure time. Within the four 
hour time window of exposure to IFN-P the transcription response 
of both 'immediate-early' (expressed independentiy of protein 
translation) and 'secondary response genes' is expected [73]. The 
assessment of the temporal nature of the IFN-P response genes was 
beyond the scope of our study; h<)we\'er, extended time course 
analyses in vitro, that may be part of future studies, could 
contribute to the identification of sustained transcription patterns 
that may better resemble the functional effect of drugs in vivo. 

Conclusions 

New IFN-P response genes were identified by foc:using on the 
response within the immune compartment of the B cells. The new 
IFN-P response genes identified by LCL analysis included genes 
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